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Abstract 

We compare the results of measurements of the magnetic susceptibility x(T), 
the linear coefficient of specific heat 7(T) = C(T)/T and 4/ occupation num- 
ber nf(T) for the intermediate valence compounds YbXCu4 (X = Ag, Cd, In, 
Mg, Tl, Zn) to the predictions of the Anderson impurity model, calculated in 
the non-crossing approximation (NCA). The crossover from the low tempera- 
ture Fermi liquid state to the high temperature local moment state is substan- 
tially slower in the compounds than predicted by the NCA; this corresponds 
to the "protracted screening" recently predicted for the Anderson Lattice. 
We present results for the dynamic susceptibility, measured through neutron 
scattering experiments, to show that the deviations between theory and ex- 
periment are not due to crystal field effects, and we present x-ray-absorption 
fine-structure (XAFS) results that show the local crystal structure around 
the X atoms is well ordered, so that the deviations probably do not arise 
from Kondo Disorder. The deviations may correlate with the background 
conduction electron density, as predicted for protracted screening. 
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I. INTRODUCTION 



The rare earth intermediate valence (IV) compoundsS are moderately heavy fermion 
compounds where the characteristic (Kondo) energy is large compared to the crystal field 
splitting (Tk > T c f). Unlike the truly heavy fermion (HF) compounds! (where Tk < T c f) 
the IV compounds do not reside close to a quantum critical point for a transition to an- 
tiferromagnetism, so non-Fermi liquid behavior is neither expected nor observed and the 
scaling behavior that is observed does not reflect proximity to such a phase transition.^ Fur- 
thermore most of the IV compounds are cubic, so that anisotropy and low-dimensionality 
are not issues. Hence the IV compounds are physical realizations of the isotropic orbitally 
degenerate (JVj = 2J + 1 = 8 for Yb) Anderson Lattice Modefl for three spatial dimensions. 
This is an archetypal problem that is both simple and elegant for exploring the physics of 
electronic correlations in solids. 

A key issue for the Anderson Lattice is the role of lattice coherence, which can be thought 
of as dispersive or band-like behavior of the 4/ electrons or alternatively as correlations 
between the 4/ electrons on different lattice sites. This contrasts the Anderson Lattice (AL) 
with the Anderson Impurity Model (AIM) where no such coherence is present. For the HF 
and IV compounds the transport behavior, which depends crucially on the periodicity of 
the scattering potential, clearly manifests lattice coherence: the low temperature resistivity 
(which is finite for the AIM) vanishes for the periodic compounds, as expected for a system 
obeying Bloch's Law;S! the de Haas- van Alphen signals are characteristic of renormalized / 
bands;! and the optical conductivity exhibits a Drude response at low temperatures that also 
reflects renormalized masses.@ On the other hand, the dynamic susceptibility x"{ u ) shows a 
Lorentzian power spectrum with very little g-dependencej suggesting that Ihe spie/valenee 
fluctuations are very local and uncorrelated, as expected for the AIM; and thermodynamic 
properties that are dominated by the spin/valence fluctuations such as the temperature 
dependence of the susceptibility x(T), specific heat C{T) and 4/ hole occupation number 
nf(T) (the Yb valence is z = 2+n/) seem to follow the predictions of the AIM, at least 
qualitatively.!'! Recent theoryS of the Anderson Lattice suggests, however, that there should 
be observable differences between the behavior of the AL and the AIM for these quantities. 
In particular, "protracted screening" can occur in the AL, which means that the crossover 
from the low temperature Fermi Liquid state to the high temperature local moment state is 
slower for the lattice case than for the impurity case. 

Protracted screening has been invokeoS to explain photoemission resultsB which show 
far less temperature dependence than expected based on the AIM. Such results are 
controversial .tJ To search for such effects in the bulk thermodynamic behavior we herein 
compare our measurements^ of x(T), C(T) and nf(T) for the series of related cubic (C15b) 
compounds YbXCu4 (X = Ag, Cd, In, Mg, Tl and Zn) to the predictions of the Anderson 
Impurity Model, calculated within a single approximation scheme, the non-crossing approx- 
imation (NCA). By making this comparison for several measurements for several related 
compounds, we put strong constraints on the applicability of the model. In addition, we 
include measurements of the dynamic susceptibility ~x!'{uS) f° r X — Ag, Mg, Tl and Zn; 
these allow us to determine whether the deviations from AIM behavior could arise from 
crystal field effects and also give a fourth experiment for comparison to the predictions of 
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the model. Finally, we include measurements of the x-ray absorption fine structure (XAFS) 
from the X-atom K edges for YbXCu 4 and X = Ag, Cd and In and the L 3 edge for X = Tl 
to determine whether deviations from AIM behavior could arise from local lattice disorder 
in the samples. 

II. EXPERIMENTAL DETAILS 

The data for the 4/ occupation number rif(T) (Fig. 1), the susceptibility x(T) (Fig. 
2) and the linear coefficient of specific heat 7 = C/T (Table 1) are taken from Sarrao 
et alM As explained in that paper, the samples were small, high quality single crystals 
grown in XCu flux. The specific heat was measured using a thermal relaxation technique. 
The susceptibility was measured using a SQUID magnetometer. In this work we have 
subtracted a small "Curie tail" (with Curie constant typically of order 10~ 2 emu-K/mol) 
from the susceptibility. The occupation number nf(T) was determined from the near-edge 
structure in Yb L3 x-ray absorption measurements. 

The neutron scattering measurements were performed in the time-of-flight mode using 
the LRMECS spectrometer at IPNS (Argonne National Laboratory). The experimental 
conditions were similar to those discussed in a recent studyEl of YbInCu4; we refer the reader 
to that publication for a more detailed account. The temperature of the measurement was 
10K, which is considerably smaller than Tk in each case. The samples utilized in these 
experiments were powders, typically of 50g mass, grown from the melt inside evacuated 
tantalum tubes for X = Ag, Mg and Zn; for X = In and Tl, a large number of single 
crystals grown as described in Sarrao et a/.0 were powdered to form the sample. The 
susceptibilities of these samples were identical to those shown in Fig. [| except for X = Ag, 
where the temperature T max of the maximum in the susceptibility was 10% smaller for the 
neutron sample, and for X = Tl where a fraction of the crystals used to create the powder 
had large Curie tails indicative of disorder or impurities. For each sample, given that the 
elastic energy resolution (FWHM) varies as 5E ~ 0.07-0.1 E iy we used several incident 
neutron energies Ei (e.g. 35, 80 and 150meV) to give greater dynamic range in energy 
transfer AE for the measurement. To improve statistics we took advantage of the lack 
of Q-dependence of the magnetic scattering^ and grouped detectors into three bins, with 
average scattering angle 20° (low Q), 60° and 100° (high Q). The scattering was put on an 
absolute scale by comparison to the measured scattering of a vanadium sample. We also 
measured the scattering of YAgCu 4 , LuMgCu 4 , YMgCu 4 , YTlCu 4 and LuZnCu 4 in order to 
help determine the nonmagnetic scattering. 

The XAFS experiments were performed on Beam Line 4-3 at the Stanford Synchrotron 
Radiation Laboratory (SSRL) using a half-tuned Si(220) double-crystal monochromator. 
Samples were ground, passed through a 30 /im sieve and brushed onto scotch tape. Pieces 
of tape were stacked such that the change in the absorption at the Yb L m edge was < 
unity. The samples were the same flux-grown crystals as used in Sarrao et alM and whose 
susceptibilities are shown in Fig. [2| For temperature control, we utilized a LHe flow cryostat. 
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III. THEORETICAL DETAILS 



The predictions of the Anderson Impurity Model were calculated in the Non-Crossing 
Approximation (NCA).@ This approximation allows for calculation of all relevant experi- 
mental quantities, both static (x, nf, 7 = C/T) and dynamic (x"{AE)), for realistic orbital 
degeneracy and for realistic spin-orbit and crystal-field splitting. It agrees wellEl over a 
broad temperature range with calculations using the Bethe Ansatz, but does produce spu- 
rious and unphysical non-Fermi Liquid artifacts in the low frequency density of states 
and in the dynamic susceptibility at very low temperatures.0 These artifacts are not sig- 
nificant for the temperatures and frequencies studied here. The conduction band was 
assumed to have a Gaussian density of states of width W and centered at the Fermi energy, 
i.e. N(e) = e~ £ l w /(y/nW). The hybridization matrix elements were assumed to be k- 
independent, i.e. 14/ = V, which may be a good approximation if the thermodynamic and 
magnetic properties are predominately determined by the conduction band states within 
/cTk of the Fermi energy. The finite temperature NCA self-consistency equations0@ were 
solved using three overlapping linear meshes. The thermodynamic and spectroscopic results 
for Ce impurities with spin-orbit splitting and in various crystal fields were compared with 
those published in Bickers et a0 and good agreement was found despite the difference in 
procedure and the different choice of conduction band density of states. 

For Yb IV compounds where crystal fields can be ignored the magnetic orbital degeneracy 
is Nj = 8 and there are essentially four input parameters for the AIM calculation: the spin- 
orbit splitting^ (which we fixed at A so = 1.3eV, the value typically observed in photoemission 
experiments^), the width W of the Gaussian conduction band, the hybridization constant 
V and the /-level energy (relative the Fermi level) Ef. Since it is only the conduction states 
within maxjfcT, kT^} of the Fermi level that contribute significantly to the experimental 
quantities discussed here, the bandwidth parameter W must be chosen to give the correct 
value of N(e-p) for the background conduction band. Assuming the background bandstates 
in YbACu4 are similar to those in LuACu4, we choose W to reproduce the linear specific 
heat coefficient 7(Lu) of the corresponding LuXCui compounds,^ using the free electron 
approximation to convert the measured value of 7(Lu) to iV(ep). This is an important 
constraint since (as we found in the earlier paper0) it is possible to force fits to nf(T) over 
the whole temperature range if unrealistically small values of the bandwidth are used in the 
model. In YbXCu4 compounds, the valence bands have a width of order 10 eVjllEil since 
the 2a full-width of the Gaussian band is ~ 4W, the values of W (~ leV, see Table U) 
obtained from 7(Lu) are quite reasonable. 

Once A so and W are fixed, V and Ef can be uniquely determined by fitting to the ground 
state values of the susceptibility and occupation number. The values of these parameters 
are given in Table |; the Kondo temperature is calculated from the formula 

V 2 1 W 6 V^WE f 

Tk = (^»' <7C>* W e ^ (1) 

which includes the effect of spin orbit splitting but ignores crystal field splitting. 

For these values of input parameters, we calculated the frequency dependence of the 
dynamic susceptibility at 10K and the temperature dependence of the static susceptibility, 
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the 4f occupation number and the free energy. From the latter we determined the linear 
coefficient of specific heat by fitting to the formula F = E — |(^-) 2 in the temperature 
range 0.03 < (T/Tk) < 0.07; the values of 7 so determined have a systematic error of order 
5-10%. 



IV. RESULTS AND ANALYSIS 

A. Neutron scattering 

The primary goal of the neutron scattering measurements was to determine whether 
crystal field excitations can be resolved in these compounds. Determination of the magnetic 
scattering in polycrystals of IV compounds requires correct subtractionof the nonmagnetic 
scattering. We have adopted a variant on the conventional procedure^ for accomplishing 
this where the YXCU4 (and/or LUXCU4) nonmagnetic counterpart to the corresponding 
YbVCu4 sample is measured to determine the factor h(AE) = S^high Q; AE;Y)/S(\ow 
Q; AE;Y) by which the high-Q scattering (which is almost totally nonmagnetic scattering) 
scales to the low Q values where the magnetic scattering is strongest. This factor varies 
smoothly from a value ~4-5, at AE = 0, to a value of 1 at large AE where Q-independent 
multiple scattering dominates the nonmagnetic scattering. (When this decrease in h(AE) 
is neglected the nonmagnetic scattering will be underestimated at large energy transfer.) In 
addition we have included an energy-independent multiplicative factor k which is needed to 
account for the fact that the multiple scattering, which dominates the low-Q nonmagnetic 
scattering, has different strengths for Yb than for Y or Lu. (When this factor is neglected 
the magnetic scattering is overestimated.) That is, we determine the magnetic scattering 
from the formula 

5 mag (low Q) = 5 tot (low Q)-k h(AE) S tot (high Q) (2) 

The magnetic scattering is related to the dynamic susceptibility through the formula 

9 N 

S mag = ^t/ 2 (Q)(1 - e^-YW, A£0 (3) 

7TfM B 

where f 2 {Q) is the Yb 4/ form factor. For each k we fit the dynamic susceptibility assuming 
that it is Q-independent and assuming a Lorentzian power spectrum: 

X »(Q,AE) = Xdc{T) AE (^){( |(Ag J n)a + n ) + ( [ (Afl + ^ + P ] )} <4) 

The prediction of the Anderson Impurity Model for x" at T << Tk also can be fit to this 
formula; we have included in Table 1 the values of E and T deduced from the fits to the 
NCA results. The fits to the experimental data include corrections for absorption and for 
instrumental resolution. The data for every Ei is included in the fit. We then find the value 
of k that gives the smallest reduced \ 2 f° r the fit. The best-fit values of k were consistently 
in the range 0.65-0.75; furthermore, the best-fit values of k gave better agreement (10-20%) 
between the fit value of Xdc and the value shown in Fig. El than was seen for k = 1. Use of 
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this multiplicative factor also leads to excellent agreement between the values of magnetic 
scattering estimated at different incident energies. In Fig. |3] we show the data (symbols) 
for the magnetic scattering at low Q (fixed average scattering angle =20°) and the fits 
(solid lines). (Plots of the data for YbInCu4 are given in the earlier publication. 0) We 
note that in the time-of-flight experiment Q varies with both and AE at fixed <fi and 
consequently the form factor, the absorption and the instrumental resolution all depend on 
Ei\ this explains why the data and the fits for different E^ don't overlap. The parameters Eq 
and T for the best fits are given in Table 1, where they are compared to the values predicted 
by the AIM calculation. We note that our results agree well with an earlier measurement 
of YbAgCu4.il Although the fits are of reasonable statistical quality with reduced x 2 — 1> 
the degree of systematic uncertainty due to the assumptions made about the nonmagnetic 
scattering is unknown. The problem is especially serious for IV compounds where Tk is 
large so that large incident energies are required, with the result that multiple scattering is 
large. To be cautious, we expect a 10% error in the determination of E and T. 



B. XAFS 

The primary goal of the XAFS measurements was to determine whether A/Cu site 
interchange is significant in these compounds. Such disorder is a real possiblity given that 
YbCus can grow in the C15b structured and given that UPdCu4, which also_grows in the 
C15b struucture, has been shown to have significant Pd/Cu site interchange.BJ Our earlier 
neutron diffraction measurement!! of flux- grown crystals of YblnCui, coupled with Rietveld 
refinement, suggest the average crystal structure is very well ordered. We employ the XAFS 
technique here both because the measurement does not depend on lattice periodicity (i.e. 
it is a local probe and is thus complimentary to diffraction) and because it is atomic- 
species specific. In such a measurement, a fine-structure function %(&;) is extracted from the 
absorption data and related to the radial bond length distribution around the absorbing A 
atomic species. This measurement is particularly sensitive to A/Cu site interchange because 
if any X atoms rest on the Cu sites, a short X-Cu bond length at about 2.55 A will appear 
that corresponds to the Cu-Cu nearest neighbor bond length in the nominal structure. This 
bond length is significantly shorter than the nearest neighbor A-Cu pairs at ~ 2.93 A or 
the A-Yb pairs at ~ 3.06 A, and so is easy to resolve if the site interchange is frequent. In 
addition to searching for A/Cu site interchange, we also look for more generic disorder in 
the nominal A-Cu and A-Yb pairs by analyzing the distribution widths extracted from the 
r-space fits as a function of temperature between 15-300 K. Bond length distribution widths 
should follow a correlated-Debye model0 as a function of temperature if the structure does 
not have any static (positional) disorder. This model differs from the usual Debye model 
because it includes the correlations in the motions between the atoms in a given pair. It has 
been shown to be accurate to within ~5% for metals.0 Any static disorder should manifest 
as a constant offset, as is clearly shown, for example, for La 1 _ 2 .Ca :r MnO3.0 

The XAFS data were fit in r-space, following procedures in Ref. [ |2~4] , |2"T|j . Initial fits 
included a short A-Cu pair at ~2.55 A. The relative amplitude of this peak compared to 
the amplitude of the nominal A-Cu at ~2.93 A gives the percentage of A atoms on Cu 
sites. Errors in this measurement were determined by assuming that the fitted statistical \ 2 
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was equal to the degrees of freedom in the fit as given by SternEI and then increasing the 
amount of site interchange until x 2 increased by one, while allowing all the other parameters 
in the fit to vary. These results are presented in Table III]. For X = In, Cd and Ag the site 
interchange has the same value as the error, and hence the results are consistent with the 
zero site interchange. For YbTlCu4 the error is smaller than the site interchange but it 
is also smaller than the error in the other cases. Furthermore, the fitted Tl-Cu bond 
length was very low (2.40 A) even though Tl has the largest atomic radius of all the X's 
measured. We note that sytematic errors in the fit can arise from use of an incorrect theoretic 
backscattering amplitudeBl Given all this, we think that the error is underestimated for 
YbTlCu4. Therefore, we conclude that, within ~3%, none of the materials measured have 
any X/Cu site interchange. 

Once the lack of site interchange was determined, we fit the data assuming no site 
interchange and determined the pair distribution widths to search for any other disorder. 
Examples of two of these fits are shown in Fig. The bond lengths of the nominal X-Cu 
and X-Yb pairs do not deviate from the average C15b structure assuming lattice parameters 
from Ref. |T^ by more than 0.02 A. The distribution widths are shown in Fig. [5], as are the 
fits to the correlated-Debye model. Parameters of the correlated-Debye fits are given in 
Table O. The fits are quite good and the observed values of the static disorder a^ tat are 
typical of values seen in well-ordered compounds. 



C. Susceptibility, 4/ occupation number and specific heat 

The results of our calculations for x(T) and n/(T) for the different YbXCu 4 compounds 
are shown in Figs. [I] and [| where they are compared to the experimental behavior. The 
most noteable feature is that for X = Cd, Mg and Zn, the experimentally determined 
occupation number n/(T) and the effective moment = Tx/Cj (where Cj is the J = 
7/2 Curie constant) rise towards the high temperature limit more slowly than predicted by 
the impurity theory. The effect is somewhat weaker but still apparent for X = Ag. For X = 
Tl the experimental occupation number shows a weak retardation compared to theory but 
the experimental susceptibility tracks the impurity theory up to room temperature. The 
linear coefficients of specific heat determined from our calculation are given in Table 1 where 
they are compared to the experimental values; the experimental and calculated values of the 
Wilson Ratio 1Z = (ix 2 R/?>Cj)x{fy/l are also compared in the table. 



V. DISCUSSION 

Figures [l] and ^| give our basic result: the crossover from the ground state Fermi Liquid 
to the high temperature local moment state is slower for the real materials than predicted 
by the Anderson Impurity Model. 

There are several reasons why this comparison to theory might be incorrect. The first 
is that crystal fields may be large in these compounds. Crystal fields can cause the ground 
state multiplet to have a lower degeneracy than the value Nj = 8 expected for J = 7/2. 
Since the AIM results depend on Nj this could lead to a discrepancy between theory and 
experiment if Nj were incorrectly assigned in the calculation. Indeed, in the earlier paperS 



7 



we found that for some cases we could fit the data assuming a lower degeneracy, e.g. J = 5/2 
for YbMgCu4. Our inelastic magnetic neutron scattering experiments were designed to test 
for this possibility. As mentioned above, there is a systematic error of order 10-20% in the 
determination of the peak position and width of the magnetic scattering. For our present 
purposes, however, the basic point is that the data can be well-fit by a single broadened 
Lorentzian whose peak energy Eq is in reasonable accord with the predicted value (i.e. it is 
of order Tk); there are no sharp crystal field levels in evidence. This also is in accord with 
the fact that the crystal fields seen in trivalent compounds such as YblnNijl or YbAuCu^Eil 
are of order 2-4meV which is significantly smaller than Tk for most of these compounds. 

A second potential problem for our analysis is that if Yb atoms are in a disordered 
environment, there can be a distribution of Kondo temperatures. In this case, the measured 
susceptibility involves an integral over this distribution of Kondo temperatures, which would 
spread the crossover over a broader range of temperature. This effect may be relevant 
in UPdCu4, where there is significant (~ 25%) Pd/Cu site exchanged Our XAFS results 
for the YbXCu4 compounds show that within ~3% there is no X/Cu site interchange; 
furthermore, the observed values of the static disorder <j 2 stat are typical of values seen in 
well-ordered compounds. Basically, the XAFS results are consistent with a very high 
degree of order for the flux-grown crystals. We thus believe that Kondo disorder is not 
a significant factor for these compounds. (A possible exception is YbZnCu4 for which we 
could not analyse the Zn if -edge XAFS because it is superimposed on the Cu if-edge XAFS. 
Given the chemical similarity of Zn and Cu, Zn/Cu site disorder seems a real possibility.) 

A third problem for this analysis concerns the interpretation of the Yb L 3 x-ray ab- 
sorption used to determine the occupation number n/(T). The analysisi'0 relies on rather 
simple assumptions: the divalent and trivalent absorption edges are assumed to have the 
same shape and the same absorption matrix element, so that n^ is determined from the frac- 
tional weights of the two features. Similar assumptions are made for other methods (XPS 
core levels; valence band UPS; Mossbauer isomer shift) for determining valence. While the 
resulting systematic error is not well-understood, it has been arguedH that the L 3 method 
is superior to use of isomer shift or lattice constant anomalies for the determination of va- 
lence. It has been pointed outE! that the values of nj estimated from the L 3 experiments are 
consistently and significantly larger than those estimated from valence band photoemission. 
The latter estimate also requires the assumption of equal absorption matirix elements for 
the di- and tri-valent configurations and is extremely sensitive to choice of background. A 
point in favor of the use of x-ray absorption is that, unlike photoemission it truly samples 
the bulk of the compound; indeed it has been arguecS that near-surface reconstruction 
gives rise to differences between the L 3 and UPS estimates of n/. Whether this is so and 
what the sytematic errors are in both measurements are, in our opinion, topics for future 
research. In any case, the deviations from the AIM that we observe for n/(T) using the 
L 3 measurement (Fig. 1) are qualitatively very similar to those observed for the effective 
moment (Fig. 2), the measurement of which is not subject to significant systematic error. 
(In our earlier papei0 we showed that at room temperature for all cases except X = Ag, rif 
is nearly equal to the effective moment.) The photoemission analysis yields an extremely 
slow crossover, to the extent that the hole occupancy is found to be nj ~ 0.6 in the high 
temperature state of YbInCu4 which is much smaller than the value (> 0.9) that we have 
observedii using L 3 . The latter value is, in our opinion, more realistic given the essentially 
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trivalent behavior of the susceptibility of YbInCu4 at high temperature. 

We therefore argue that we have ruled out most of the major sources of error in our 
analysis, so that we can have a reasonable degree of confidence that the slow crossover is a 
real effect. As mentioned in the introduction, recent theory0 for the (S = 1/2) Anderson 
Lattice, treated in the 1/d approximation (where d is the spatial dimension of the system) 
predicts such a slow crossover, referring to it as "protracted screening". The degree of 
slowness is sensitive to the filling of the background conduction band; the crossover is slowest 
for partial filling of the background band. At present it is unclear the extent to which the 
effect depends on the orbital degeneracy. 

The deviations of the data (Figs. 1 and 2) for YbACu4 from the predictions of the AIM 
are strongest for X = Mg, Zn and Cd, weaker for X = Ag and weakest for X = Tl. We note 
that there is no correlation between this trend and the magnitude of the Kondo temperature 
(Table 1) nor of the ground state occupation number (Fig. 1). There does appear to be 
a correlation with the magnitude of the Hall coefficient for LuXCu^ Rq = —1.8,-1.0 
and -0.8 x 10- 10 m 3 /C for X = Mg, Zn, Cd;@ -0.6 x 10- W m 3 /C for X = Ag;0 and 
< 1 x 10- n m 3 /C for X = T10 This quantity can be taken as a measure of the conduction 
electron density in the background band. In a one-band model, where Ro = 1/ne and 
where there are 24 atoms in a unit cell of side 7.1 A these values of Ro correspond to the 
values 0.52, 0.93, 1.16 and 1.55 electrons per atom for the sequence {X = Mg, Zn, Cd, 
Ag}. While the one-band assumption is unrealistic, nevertheless these results suggest that 
the deviations from the AIM become weaker as the carrier density becomes larger. This 
accords with theory0 which predicts no deviation for half filling of the background band 
and increasing deviation as the conduction electron density decreases. 

We should point out that a slow crossover is expected on rather general grounds. In 
the single impurity calculations the position of the chemical potential /x is set by the host 
metal conduction band density of states, and the chemical potential is independent of tem- 
perature to order kT/W. In the hypothetical case of a periodic array of impurities that 
are noninteracting and incoherent the appreciable /-weight and its asymmetric temperature 
dependence will lead to a significant temperature dependence of /x(T) and hence lead to a 
deviation of nf(T) from the behavior found in the solution to the single impurity model. 
It is not clear to us the extent to which this effect contributes to the predicted protracted 
screening. 

In addition to our main conclusion, we would like to point to two more features of the 
analysis. The first is that the experimentally determined linear coefficients 7 of specific 
heat for X = Mg, Cd, Ag and Zn are 20-30% larger than the values predicted by the AIM. 
This means that the experimental value of the Wilson ratio 1Z is smaller than predicted 
by the AIM. This has been noticed before. !0 For YbInCu4 there is excellent agreement 
between the experimental and calculated values; for YbTlCu4 the calculated value of 7 is 
30% smaller than the experimental value. The second feature is that there is considerable 
disagreement between the measured and predicted values of the neutron lineshape parame- 
ters Eq and T given in Table 1. We remind the reader that the degree of systematic error 
for this determination is poorly understood, due to the difficulty of determining the non- 
magnetic scattering for large Tk compounds. In addition, all Q dependence of the lineshape 
in polycrystals is automatically averaged out. In the case of YbInCu4, we have confidence 
in the results, since an alternate means of subtracting the nonmagnetic scattering^ agrees 
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with the method used here and since a neutron scattering study of a single crystalQ showed 
that the Q-dependence of the lineshape is quite small. For that case we note that the ex- 
perimental peak energy E agrees well with the prediction of the AIM but the measured 
linewidth T is a factor of two smaller than predicted, as though the spin fluctations have 
a longer lifetime than predicted by the AIM. On the other hand for YbMgCu4, YbTlCu4 
and YbZnCu4 the experimental values of peak energy differ significantly from the calculated 
values and the linewidths are larger than the predicted values. Such a large linewidth could 
be a consequence of a Q-dependence of the peak parameter E (Q); studies of single crystals 
are necessary to establish whether such a Q-dependence is present. Finally we note that if 
the lineshape parameters for YbAgCu4 are increased by 10% to account for the 10% smaller 
value of T max observed for the neutron sample, they are in excellent agreement with the 
AIM predictions. 

In conclusion, we have given evidence that the crossover from low temperature Fermi 
Liquid behavior to high temperature local moment behavior is slower for periodic intermedi- 
ate valence compounds than predicted by the Anderson Impurity Model. We have included 
supporting experiments to rule out the possibility that the disagreements between the data 
and the AIM predictions are due to large crystal fields or to Kondo disorder. We feel 
that the results are fairly robust. It will be interesting to see whether similar "protracted 
screening" is observed in other compounds, especially in those based on Ce where the orbital 
degeneracy is smaller (Nj = 6) and where the size of the 4/ orbital is considerably larger. 
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TABLES 



TABLE I. Input parameters W, Ef and V for the AIM calculation; the calculated Kondo 
temperature Tk; and the theoretical and experimental values of the specific heat coefficient 7, the 
Wilson ratio 1Z and the neutron lineshape parameters Eq and L. 



Compound 


W 
(eV) 


E f 
(eV) 


V 

(eV) 


T K 
(K) 


7 

( TflJ \ 

\mol-K 2 > 


TZ 


Eq 

(meV) 


L 

(meV) 


YbTlCu 4 


1.286 


-0.50146 


0.2195 


514 


32.1 


1.26 


54.0 


31.2 










Expt: 


24.2 


1.67 


38.8 


32.3 


YbMgCu 4 


1.005 


-0.1897 


0.128 


500 


36.7 


1.22 


44.6 


23.4 










Expt: 


53.3 


0.84 


31.9 


34.5 


YbInCu 4 


l 1 


-0.7442 


0.232 


299 


41.3 


1.44 


41.6 


25.9 










Expt: 


41.3 1 


1.44 1 


41.6 


13.5 


YbCdCu 4 


0.929 


-0.17576 


0.098 


127 


120.4 


1.40 














Expt: 


165.6 


1.02 






YbAgCu 4 


0.865 


-0.4485 


0.148 


95 


137.9 


1.29 


11.3 


7.1 










Expt:: 


198.9 


0.89 


9.8 2 


6.0 2 


YbZnCu 4 


1.213 


-0.1984 


0.1038 


60 


296.8 


1.33 


5.0 


4.1 










Expt: 


370 


1.07 


7.4 


6.1 



Notes: 1) The values of W are deduced from the specific heat of the corresponding LuXCu4 
compound and the experimental specific heat coefficients have been corrected for the non-4/ 
contribution by subtracting the specific heat coefficient of the LuXCu 4 compound. Since 
the groundstate of LuInCu4 is semimetallidil and does not serve for these estimates, we have 
arbitrarily set W = leV and we have subtracted a corresponding amount (8.7mJ/molK 2 ) 
from the specific heat coefficient. 2) For comparison to theory these values should be 
increased by 10% to account for the fact that the neutron scattering sample had a Kondo 
temperature that is 10% smaller than the sample used to measure x an d fif. 



TABLE II. XAFS fit results between 2.0 and 3.2 A. Fits only include the nominal structural 
peaks X-Cu at ~2.9 A and X-Yh at ~3.06 A. S$ is an amplitude reduction fact that sets the scale 
for the distribution variance measurements. See text for details. 







e cD (K) 


2 

^stat 


(A 2 ) 


X/Cu 


X 


^0 


Cu 


Yb 


Cu 


Yb 


interchange 


Tl 


0.89(5) 


230(5) 


230(5) 


0.0005(4) 


0.0005(5) 


4(1)% 


In 


1.04(5) 


252(5) 


280(5) 


0.0009(4) 


0.0011(5) 


2(3)% 


Cd 


0.98(5) 


240(5) 


255(5) 


0.0007(4) 


0.0010(5) 


5(5)% 


Ag 


0.91(5) 


250(5) 


235(5) 


0.0008(4) 


0.0006(5) 


2(2)% 



13 



FIGURES 



FIG. 1. The 4/ hole occupation number n/(T) versus temperature for the YbYCu4 com- 
pounds. The open circles are the experimental data and the solid lines are the predictions of the 
Anderson Impurity Model (AIM), with input parameters given in Table 1. 

FIG. 2. The magnetic susceptibility x(T) (open circles) and the effective moment 
/Xgjj = Tx/Cj (solid circles) for the YbYCu4 compounds. The solid (dotted) lines are the predic- 
tion of the AIM for the susceptibility (effective moment). 

FIG. 3. The inelastic magnetic neutron scattering spectra for the YbYCu4 compounds at T 
= 10 K. For each sample, two or three values of the incident energy Ei were used to increase the 
dynamic range. (YbMgCu4: solid circles, Ei = 150 meV; open circles, Ei = 80meV; open diamonds, 
Ei = 35meV. YbTlCu4: solid circles, Ei = 150meV; open circles, Ei = 80meV. YbZnCu4: solid 
circles, Ei = 80meV; open circles, Ei = 25meV; open diamonds, Ei = 15meV. YbAgCu4: solid 
circles, Ei = 60meV; open circles, Ei = 25meV.) The solid lines represent fits to the Lorentzian 
power spectrum (see text) with the parameters given in Table 1. 

FIG. 4. XAFS data and fits on (a) YbAgCu 4 and (b) YbCdCu 4 . The outer envelope shows ± 
the amplitude (or modulus) and the oscillating inner curve is the real part of the complex Fourier 
transform of k 3 x(k)- The r-axis includes different phase shifts for each coordination shell, so that 
the r of a given peak appears somewhat longer than the actual pair distance R. These transforms 
are from 2.5-15 A 1 and Gaussian narrowed by 0.3 A -1 . Fits are from 2.0 to 3.2 A. 

FIG. 5. The variance in (a) the Y-Cu pair distribution at ~2.9 A, and (b) the Y-Yb pair 
distribution at ~ 3.06 A. The temperature dependence follows a correlated-Debye model well for 
all the data, as shown. Most of the low temperature data are very similar, except the Tl-Cu and 
Tl-Yb cr 2 's are somwhat lower due to the greater mass of Tl compared to the other Y's. 
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